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INTRODUCTION 


In this report, two subjects are investigated: 

I. Analysis of radiation damage of MOSFET data from 
Explorer XXXIV (IMP-F), and 

II. Radiation damage characteristics of MOSFET with boron 
diffused between silicon semiconductor and silicon oxide. 

The first subject is an interpretation of the discrepancy 
between the space data and the laboratory data. Accroding to 
the work of John L. Wolfgang, Jr. 1 ^ for an equal flurence, the 
space damage is less than the damage observed in the laboratory. 
Present work will establish a mathematical model based on which 
to calculate the accumulated damage especially when the radiation 
flux is not a constant, such as in the space orbit: a satellite 

would see the flux as a function of time. This result will be 
useful in the forecasting of the stability of MOSFET in 
radiation environment. 

The second subject is an attempt to analyze the radiation 
damage characteristic of MOSFET when there is modification of 
the electrical property in gate oxide region. In such a system. 
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the recent result of Vitaly Danchenko*' indicates a spontaneous 
recovery of the radiation damage. We will propose a mechanism 
for this recovery. A substantiation of a mechanism which leads 
to a spontaneous recovery in MOSFET would be very important to 
the design of MOSFET for future space applications. 

(I) Radiation Damge of MOSFET in Space 
We quote a summary of Wolfgang" 

"After 60 orbits of IMP-F 9.5 x 10 1U electrons/cm 
(energy 0.55 Mev) dose has been received by the 0.25 
gm/cm 2 shielded devices. The gate threshold shift of 
similar devices, at the same dose of 1.5 Mev electrons 
in laboratory studies, indicating possible minor anneal- 
ing is occurring over the highly elliptical orbit 
in flight." 

There are two problems in order to discuss this interesting 
observation: (1) can there be annealing at the temperature 

range of 10 to 20°C where the device experienced in IMP-F. The 
answer is affirmative, at least in the qualitative sense. A 
small annealing at a low temperature of 25°C is reported by 
Danchenko et al.^ with an observational period of 100 hours. 

This time range is close to the orbital period of IMP-F. 
Independently, we will also prove that the space data implies 
the existence of an annealing of MOSFET in space. 

In the experiment of Wolfgang in the laboratory, the 
radiation time is in the order of minutes. To reach the same 
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level of-. fluence in space will take several months. Therefore, 
our second problem is to investigate theoretically, by intro- 
ducing explicitly the annealing effect, the radiation damage 
under these two condtions. 

In the following, we will establish a kinetic equation 
for simultaneous radiation and annealing, and from this equation, 
we will discuss the observational results. 

A kinetic equation describes the change of the defect 
concentrations of a system. In low level, moderately high 
energy radiation, such is the case under present investigation 
n is linearly proportional to the fluence $. Furthermore, we 
reploted the data of Wolfgang on the fluence dependence of 
the shift of gate threshold voltage AV^, instead of seim-log, 
in a regular graph, the result shows a straight line (Fig. 1), 

f 

that is AV gt is linearly proportional to $>• Therefore, 

Vqij, is linearly proportional to n. 

V GT = 2.1 x 10“ 9 mv/(e/cm 2 ). In other words, (1) 

AV gt is related to n by a constant multiplier k, i.e., 

AV GT = kn (1') 

In our treatment, therefore, we can simply discuss the 
conventional variable n of the kinetic equation, instead of 
the complex quantity, A ^QT’ 

If n Q denotes the accumulated concentration of defects, 
in our case, as we have already discussed, n Q is linearly 
proportional to the total fluence of the radiation, 

n Q = a (proportionately constant )•$ (the fluence). (2) 
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Now, if n anneals according to the usual first 
order/ kinetics, at time t after the radiation, 
n(t) = n Q e~ t/r 

where t is the recovery (decay of defects) rate. 

In the space, since the radiation and the recovery are 
simultaneous, therefore, 

St - + 

where <J> is the flux rate. For constant flux rate, 

4 > = Therefore, starting with zero damage, after a 

time t, 

n = ctT(J>(l-e ^ 

- "o-T (l-P‘ tA) 

The result of Eq. (3) is obviously different from that 
of Eq. (5) unless in the limiting case where t approaches 
zero. The result of Eq. (4) is based on an assumption that 
I is a constant, independent of the time. We now treat the 
case that J is a periodic function of time and to be concrete 
in the form 


(5) 


4> = Asin wt 


( 6 ) 


This form of flux function has been applied by Dr. R. Waddel 
in his analysis of solar cell radiation damage for a series 

4 ) 

of Relays and ATS satellites and found a good agreement. 

In the case of Explorer, the spacecraft is in the radiation 
belt for only 4 hours orbit, i.e. about ^ of an orbit. 

Eq. (6) will represent a better approximation than the 
constant flux treatment of Eq. (5)* From Eq. (6), after 
one period, we can calculate A of Eq. (6) by a comparison 
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with the fluence $>: 

$ 



Asln tot 


Air 

2(o 


( 7 ) 


The integration of Eq. (4), now becomes, 

n = ce" t//r + ct$” cos2iot + 2(OTSin2(ot . ](g) 

1 + 4 co 2 t 2 


where a is the defect production rate. Since the unit 
of integration constant c in Eq.(8) and a are unspecified 
so far, according to Eq. (1 ! ), n of Eq.(8) can be replaced 
by AV^ which, for simplicity, will be denoted by V. 

At t = o, i.e., before irradiation, V = o, from which 
the integration constant c can be determined and the integ- 
ration from o to ir/to gives, 

V 1 = SLf- (l-x)(l-y), (9) 

where the abbreviation symbols are, s = tt/l ot, 
x = (1 + 4 co 2 t 2 ) -1 , y = e~ s . Nov;, in the laboratory 
experiment, radiation time is negligible, therefore 
Eq. (2) becomes. 


= a$y, 

where the superscript ( ' ) of V denotes the laboratory 
radiation value of V. From Eqs . (9) and (10), 

t 

V-, 

= s(l + — )/(e s -l) 

V 1 


( 10 ) 


(ID 


By a series expansion, the righthand side of Eq. (11) 


«v 
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becomes , 


(i + -si)/(i + f + {£ + ii +...), 

4tt 2 2 b 2 


The coefficient of s 2 in numerator^ (jpjr) » is smaller than 

that of denominator ( 7 ) , and all other terms in the 

*6 

denominator is positive, therefore. 


> 1. 


( 12 ) 


This result is opposite to the observation of Wolfgang. 

The above investigation is based on the result of 
a single orbit. A practical problem is to investigate the 
resultant damage after a satellite is passed through the 
radiation belt many times. 

Prom Eq. ( 8 ), for the second orbit, at t = o, V will 
have the value V.^ given by Eq. (9). At the end of second 
period, we find, from the same model of calculation, 

V, - (1-xKl-y 2 ). (13) 

and at the end of k-th period, with a method of deduction. 


V^. = “ (l-x)(l-y k ) 


(14) 


On the otherhand, in the laboratory experiment, 
= f(y 2 +y) and 


V k * 


fy (l-y k )/(l-y) . 


(15) 


The ratio between and is. 


1— x 1-y 

r = — — 

s y 


(16) 
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That is ? r is independent of k. The Explorer result shows that 
r becomes smaller as k becomes larger. Therefore, we proved 
conclusively the first order kinetic model can neither des- 
cribe the radiation damage of MOSPET in a single orbit, nor 
in an accumulation of orbits. 

A similar situation has been recognized in a different 
electronic device. In the radiation damage and annealing of 
s'olar sells, we found a fractionation effect, namely, the 
annealable damage is dependent on the instantaneous degree of 
the damage, and if radiation is interrupted with an intermittent 
annealing, the resultant damage is much less than the damage 
caused by a continuous radiation with the same flurence, over 
a time period including both that of radiation and of annealing. 
We found that the result can be explained by assuming that 
the activation energy of annealing increases as the degree of 
the damage is increased. ^ The modified kinetic equation 
thus becomes, 

dn = _ exp(-6n 1/3 ) + (17) 

dt t 

0 

1/3 

where x of Eq. (4) is now replaced by x Q exp( n J ), 

with x the recovery time at infinitesimal damage, and 3 
o 

is a numerical parameter related to the radiation sensitivity 
of the device . 

To investigate Eq. (17) analytically is not possible 
because of its highly non-linear nature. Therefore, numerical 
computation has to be applied and this problem will be presently 


discussed. 
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In Eq.(6), if t is measured in the unit of orbital 
period, then a = tt Denoting t = t o ii, and g 3 n = m, Eq. (17) 
is reduced to 

- mexp(-m 1/3 ; - e 3 r o aASin 2 (irT o vi) (18) 

that is, the kinetic equation has two independent parameters, 

A and ttt o . In addition, since we are only interested in 
the long term behaviour, r x q p will represent many orbits and 
the detailed behavior of each orbit is not an immediate 
interest. Therefore, Eq. (18) is only a one-parameter equation, 
namely, that of |3 3 t q uA which will be denoted by A*. Since 
m is a linear function of n, m is a linear function of Vq^V, 
one can directly study V as a function of B. 

The method to determine B from experimental data is 
as follows: 

(1) With the best estimation from the data, V shows 

a periodicity with the same period as that of the satellites, 
and the absolute amplitude of the periodic V, that is, from 
the minimum to maximum of AV^ at each period is about 
20 mv, the quantization level of analog-digital conversion. 

(2) The slope of V is a function of t, i.e. the number 
of orbit, after a long time in orbit is almost a straight line 
with a slope : 1 mv/orbit. 
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The above quantitative data are for a special MOSFET 
designated as VGT1 by Wolfgang. Different transitors 
behave in a quantitatively different way. This can be 
accommodated by changing the values of parameter A'. Since 
A is the same for all transistors as they all experience same 
radiation environment, the cause of variation in A* among 
different transitors is caused by 8,t q or a and the result is 
very sensitive: a small change in A' causes a drastic differ- 

ence in the behaviour of m. Fig. 2 illustrates this point. 
When A* = 2.5, there is a rapid increase of m in the early 
orbits, and then reaches a saturation. When A’ value is 
doubled, m increases almost linearly as the number of orbits 
is increased. 

Fig. 3 shows in detail the function m for three cases, 
all under the same sine square radiation orbit and a 
periodic fluctuation is observed in all cases. Curve (1) is 
for first order kinetics. As can be seen from Eq. (8), m 
saturates at large number of orbits, independent of the value 
of A*. If t is infinitely large, that is, there is no 
annealing term, .increases with periodic inflections, but 
without minima, and divergences. In the case of our proposed 
Eq. (18), m has initial rapid rise, but with obvious maxima 
and minima, eventually reaches a saturation or a continuous 
rise, as discussed in Fig. 2. 

We will show now that Eq. (18) represents the space 
data well. First, we will discuss the data of the transistor 
VGT8. Wolfgang fit the data empirically by a straight line 
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without physical interpretation. The curve obtained from 
Eq. (18) shows a better fit in the sense of smaller standard 
deviation. The value of AV^ is obtained from a comparison 
of AV QT and m at 140 orbits which gives m=1.0AVQ,p(Fig 4.) 

The next data we wish to discuss is for VGT1. By fitting 
at the data near 140 orbits, we obtained the conversion con- 
stant from m to for different values of A'. (See Fig. 5) 

When A* changes from 0 to 2.5, a satisfactory fit is 
obtained when A' = 2.8. A remarkable result is an initial 
rapid increase of AV Gr p. This effect was observed in space 
data and Wolfgang thought this was due to a large temperature 
variation. A difficulty in this interpretation is that 
taansistors with the same circuit configuration show a 
same behavior with good statistics, and is different appre- 
ciably from other transistors with different circuitry con- 
figurations. Therefore, it is not possible to explain the 
observed permanent effect. We have been able to show here 
this large shift is an inherent feature of the kinetic 
equation and thus, strengthened the validity of our proposed 
annealing model. 

(2) Diffusion Theoretical Approach to Annealing Due to 

Acceptor Impurity Doping 

I. Introduction 

Radiation produces an interfacial positive charge state 
in metal-oxide semiconductor system and resulted in a shift 
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of the gate threshold voltage. The positive charge can 
be annealed if it diffuses toward the semiconductor region 
where the charge is more mobile and the medium is conductive, 
therefore, there is no electric field effect. This process 
could be accelerated in principle ; 'by an impurity drift field 
with proper sign of charges and the charge distribution. 
HOwever, an intentional attempt to utilize this effect has 

O \ 

been carried out only in a limited way. In particular, 
there is some difficulty to duplicate the specimen, and we 
Intend to conclude that the stringency of the charge distri- 
bution may be beyond the control capability of producing the 
desired devices. In this respect, a theoretical analysis based 
on the mathematical model which reappear some physical model 
could be rewarding V ■ in the economy of experimental effort 
and the understanding of the problem. 


II. Diffusion Model 

In the previous report, we outlined one such approach 
based on diffusion equation 'Of following type: 


- VE 


3n 

3 X 


n 


( 1 ) 


where n is the density of radiation induced positive charge. 

D is the diffusion coefficient, y is the mobility, and t is 
the life time. E is the field due to a charge distribution: 
we have assumed an exponential concentration profile and 
therefore E is independent of x. We have solved Eq. (1) by 
standard procedure in previous report and the result gives, 


C 
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t ^ \ 2 n 0 , gfioo - fl rn + 

n(o,t } = T|2 e ' zq| s 2 (B+k) e m {[(B+k) + 1]. 

m=-°° 

( 6+k) 6 + 0 2 }cose (2) 

m m 

where k is the surface recombination constant, 

TT 0 

3= — with z the diffusion length, and 6 = d/Jl, d is 
kT 

the thickness of the interface region. The charge density 

at the termination of radiation is n , that is, we treat 

o 

the radiation and recovery as separate steps. Two most 

important quantities are, 9 and n : 9 are the roots 

m m m 

of the equation, 

e m + 6(6+k ) = °. (3) 

Prom the nature of trigonometrical function, m are 
integers extending from - » to + ®. 
xri m = 1 + 3 2 + 

From Eq. (4), if 6 is sufficiently small, the exponential 
series converges rapidly when t is large. When this 
very small, the summation becomes very cumbersome. Pre- 
viously, we solved small t case by a different series expan- 
sion which involves successive derivatives of transcendental 
functions. This turns out to be quite impractical for 
computer calculation and this problem is yet to be solved. 

Now Eq. (2) is solved for a ’'pulsed" sources. In the 
present configuration, the defects are annealed with an 
accumulative effect and therefore, from Eq. (2), 
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n(t) 


n. 


( n(o,t') dt’ 


n 


2n 


o —8 6 
— e 


I 

m =-0» 


0 


m (g+k)(l-e 


njnt )/{C(B+k)6+ 1]. 


(g+k) 6 + e*} n m cos0 


m 


Numerical calculations have been made on different set 
of values of g, k and t with the time in the unit of x . . 

Pig. 6 shows one example where the parameters are, 

2kT 5kT 

= 0.193 at 300°k and V = 10 mv, 
k = .1, 6=2. 

This is for Curve I in the figure . Curve II is for same 
values of k and 6, but 8=0. The total annealed defect in 
the presence of the drift field (curve I) is obviously higher 
then the case of no drift field. 

The above results are obtained with one assumption, 

namely, the concentration of defect is very large such that 

the defects removed by annealing is negligiable. The case 

for finite concentration is not carried out in the present 

work. The above results also assume that the diffusion of 

defects through the interface is an uncorrelated diffusion. 

The formulation for correlated diffusion can be found in 

Q ) 

the work of Waite However, in this approach, one can only 
introduce impurities xdth uniform distribution, that is. 
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there Is no field effect. Similar results apply to the 
model of Simons and Hughes' 1 " ' where the spatial dependence 
such as barrier effect are not possible to be introduced and 
one wonders under such a restriction, whether the physical 
model has suffered. 


SUMMARY 

This report covers two subjects and the principal 
results are: 

1) For the analysis of radiation damage data of MOSFETS, 
we proposed a kinetic equation for simultaneous radiation 
damage and annealing. The equation fits well the data from 
Explorer 34 reported by Wolfgang. Therefore, besides providing 
physical insight to the behaviour of MOSFETS in radiation 
environment, the results provide a theoretical basis to predict 
the long term behaviour of MOSFET devices in space satellites. 

2) For a modified MOS structure by introducing boron 
layer, we formulated a diffusion model which supports the 
idea that a concentration gradiant field can be introduced 
which would accelerate the annealing of radiation damage. 

I would like to acknowledge the kind cooperation of V. 
Danchenko and J. Wolfgang in discussing their interesting 
experimental results and I sincerely hope that this technical 
cooperation would continue beyond the termination of the con- 
tractual relation. Numerical calculations on computers were 
made by E. Richards and I am indebted to him for the results. 
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LIST OP FIGURES 

1) Laboratory calibrated shift of gate threshold voltage 
as a function of the fluence of (Mev electrons (reported 
from data of J. Wolfgang) 

2) Solution of Eq. ( 1 8 ) for two values of A’, the bars along the 
curve indicate the amplitude of periodic fluctuations. 

3. Solution of three kinetic equations. Curve (1) is for Eq. (4) 
- Curve (3) is for Eq. (18); Curve (2) is for Eq. (18) with 
t », i.e. zero annealing. All are evaluated for 

Eq. (6) with A’ = 1. 

4) The shift of gate threshold voltage of VGT6 : data are 
circles, broken line is an empirical fitting by Wolfgang, 
solid line is from Eq. (18) with A’ = 3 and m = 1.0 x 

a ^GT* t o = 0.5. The bar along the solid curve indicates 
the .amplitude of periodic fluctuations. 

5) The shift of gate threshold voltage of VGT1: data are 

circles. Curve (1) is an empirical straight line fitting. 
Curve (2) is for A’ = 3 S and AV^ - 3.5 m; Curve (3) 

is for A ! - 2.8 and AV^ = 8.0m. Curve (4) is for A ? 

= 2.5 and AV^ = 17.2 m. All t q = 0.5. 

6) Time rate of disappearance of radiation induced defect, 

(I) with a boron gradient field, (II), without the field. 
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